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Due to the growth of electric end-uses, the management of the variations in time of the electric power 
demand has become essential, especially in the residential sector. According to this issue, the 
anticipation of the power demand is of great interest. This implies a better knowledge of the electric 
load curve of the household stock. Papers about understanding and forecasting energy demand are 
numerous but studies on building’s load curves are rare. In this paper we propose a cross analysis of 


Keywords: 

Load curve modelling 
Literature survey 
Classification 
Residential sector 


some existing methods capable of building up a residential electric load curve. Two main types of load 
curve models have been identified in the literature: top-down and bottom-up methods. Even if the 
review presents two existing top-down approaches, the authors focused the further analysis on 
bottom-up models. For each of them we first identify its functional characteristics: finality and scope, 
input data required, output format, modeled appliances and end-uses covered, generation of the 
diversity and validation of the model. Secondly, we establish a bloc diagram representing its 
architecture with focus on the mathematical model chosen. Finally, the authors list the limits of the 
model in view of the criteria needed to build up an ideal, bottom-up and technically explicit load curve 
model for the residential sector. 

A cross reading of the different methodologies is proposed through a global table that characterizes 
and sums up the analyzed models. Moreover, a graphical representation of the models studied is 
proposed, according to three criteria (range of application, modelling of the diversity, time scale 
accuracy) that allows us to compare them “at a glance”. To conclude, an identification of the gaps 
among existing models is proposed. It consists of listing the end-uses, the appliances and other 
behavioural correlations which may affect the load curve and that are not included in the reviewed 
methodologies and need further research. 

© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Modelling the electricity consumption of the residential sector 
is not a recent research topic. For a few decades authors have 
published their methods developed to anticipate the temporal 
evolution of the electric energy demand with different time and 
space scale considerations. Among the motivations for developing 
such a model, the main one is the need of a tool to better predict, 
quantify and plan the future requirements in terms of power plants. 
Other reasons are the electric network control and management or 
the impacts of various evolution scenarii that focus on technical, 
demographic, behavioural or economical aspect. 

For such purposes the authors generally focused on the annual 
electric energy demand and the chosen methodologies are differ- 
ent depending on the finality, the scope, the input data and the 
output format of each model. 

Historically, the methods have been divided into two cate- 
gories: top-down and bottom-up approaches. Swan and Ugursal 
indicate that this classification “is with reference to the hierarchal 
position of data inputs as compared to the housing sector as a 
whole” [1]. Models operate according to one of the processes 
represented in Fig. 1. However, some authors tried to benefit from 
the advantages of both approaches and thus develop hybrid 
models. 

Swan and Ugursal [1] have classified the domestic electricity 
demand models of the literature with regard to the method used. 
They established a finer segmentation of the available methods 
than the simple differentiation between top-down and bottom-up 
methods. In this paper, we carried out a similar, global and 
transverse analysis of the methodologies for modelling the 
electric power demand. 


Bottom-up models Top-down models 


fooooo] macroscopic data 


Housing 
| stock 


Housing |€ 


= T ai | stock oo 
= =| ajaja 


microscopic data 


calculation of an energy consumption 

for a dwelling or a group of dwellings 

and extrapolation to the total housing 
stock 


attribution of an energy consumption 
to the housing stock according 
to its characteristics 


Fig. 1. Top-down vs bottom-up models. 


The power demand issue is a more recent research topic than 
the study of electric energy consumption. However, many refer- 
ences exploring the load curve research field have been published 
since the 1940s: Hamilton [2,3] and Bary [4] respectively focused 
on the properties of aggregated load curves and on the mathe- 
matical relationships between load curve characteristic factors 
(coincidence factor, load factor). In Europe, UNIPEDE [5] pub- 
lished the first load curve manual to clarify the load curve 
concepts and terminology. Rehberg [6] proposed a publication 
about the study of the residential load curve. In the United States, 
the Association of Edison Illuminating Companies (AEIC) [7-9] 
published several versions of the procedures to carry out load 
curve surveys. 

In fact power demand is more complicated to predict than 
energy demand because of its random nature and its acute 
fluctuating aspect. On top of that, even if customers can have 
the same general characteristics (e.g. two domestic customers 
with the same household size, type of dwelling, building char- 
acteristics, appliances ownership, etc.) there is every chance that 
their corresponding power demand, so their load curves for 
selected days, are completely different. Diversity that represents 
both the non-coincidence in energy use and an unlimited variety 
of customers’ characteristics is responsible for this effect. 

The influence of the human behaviour on the domestic power 
demand is so important that there is every chance for instance 
that two households with the same daily energy consumption 
will not show a similar load curve. That is the reason why 
modelling the power demand in a specific area for a specified 
target is also an hard task. 

The power demand issue was explored to ensure the safety 
of electricity supply, to better predict the peak power demand, 
to optimize the network control techniques and to analyze 
the impact of Demand Side Management (DSM) strategies or 
the modification of the network load flows after integration of 
renewable energy sources. 

The aim of this paper was the establishment of a new and 
improved classification of the residential load curve models 
described in the recent publications. This segmentation can be 
seen in Fig. 2. 

We decided to analyze these models with reference to an ideal 
and to our knowledge of non-existing model. To our opinion such 
an ideal model should present the following characteristics: 


e it has to be parametric in order to simulate various scenarii; 

e it has to be technically explicit i. e. the different specificities 
of the simulated elements (equipments, buildings, etc.) must 
specifically impact the load curve calculations and results; 
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Fig. 2. Proposed classification of the load curve models. 


e it has to be evolutive, i. e. new elements can be introduced in 
the model so as to be simulated; 

e it has to be aggregative so that results can be obtained at 
different levels (household, area, city, region, etc.); 

e all domestic end-uses must be considered in the load curve 
calculations: Heating, Ventilation and Air Conditionning 
(HVAC), Domestic Hot Water (DHW) and specific electricity 
appliances. 


According to the manner diversity is gotten, we define five 
types of models: 


e deterministic statistical disaggregation model: this approach 
consists of disaggregating measured load profiles to identify 
various appliances. Diversity is not modelled at all because it is 
embedded, in a deterministic way, in the measured data; 

e statistical random model: in such a model, the reconstitution 
of the diversity makes use of statistical data. To generate 
variations for a given scenario, a random procedure is 
applicated; 

e probabilistic empirical model: from real collected data con- 
cerning, among other things, domestic habits of people, 
probabilistic procedures are defined and applied to generate 
a diversity of results; 

e Time of use (TOU) based model: diversity is constructed 
thanks to data coming from time of use surveys (real and 
precise data concerning the behaviour of people); 

e statistical engineering model: in this case, the diversity is 
partly embedded in the measured data that serve as input 
(dwelling characteristics, weather data, penetration rates, etc.). 
On top of that, diversity is embedded in the statistical 
coefficients that adjust the original results. These coefficients 
are calculated with the help of measured load curves and 
socio-economic data. 


To carry out our analysis, we decided to follow a generic 
reading process for each model. We will focus on a series of 
criteria in order to build a standardized ID for each approach. This 
will help us to set the cross reading table at the end of this paper, 
all models considered. 

In Section 2, after a short introduction about the top-down 
approach, the authors review two of them that use a deterministic 
statistical disaggregation method. Section 3 first provides a short 


description of the bottom-up methodology, then an analysis of 
the corresponding models. These methodologies are classified 
according to the proposed segmentation. In the following, Section 
4 exposes a model which is not sortable into the two main 
families: it is an hybrid method. At the end of this paper, a 
comparison of all methods, following a cross analysis procedure, 
is carried out. The main properties of the models are put together 
into a synthetic table and a graphical representation using three 
criteria is drawn. 


2. Top-down models 
2.1, Introduction 


Top-down methodologies consider situation on a whole (making 
use of national energy statistics for instance) and try to attribute an 
electricity consumption to the studied household stock with regard 
to its characteristics. The input data for this type of model are very 
general information such as: 


e Gross Domestic Product (GDP), 

e unemployment rate, 

e present statistics on the targeted population with possibly 
predicted evolution, 

e appliance saturation rates, etc. 


Generally these methodologies consist of a mathematical 
identification of the past and present electricity demand for a 
specified area. On this basis, the model may be used to generate 
trend evolutions: that can be well predicted by this type of model. 
The advantages of such a method is its simplicity, the use of 
generally well available and aggregated data and the possibility to 
calibrate the simulation tool with real electricity consumption 
values. Yet these models are unable to simulate the impacts of 
technical breaks or the appearance of new end-uses which have 
not been measured. 


2.2. Deterministic statistical disaggregation models 


In this section we present two methodologies whose aim is the 
analysis of residential load curve. In fact, total load curve (i. e. the 
aggregated power demand of a dwelling, all end-uses being 
considered) measurements are used as input data and the 
models try then to differentiate the individual contribution of 
the observed end-uses. This work corresponds to the load curve 
analysis and it makes use of mathematical techniques such as 
regression, statistical or econometric methods. In this approach, 
the authors try to explain the domestic electric power demand 
(dependent variable) with the help of some various and indepen- 
dent variables. 


2.2.1. Aigner et al.’s model 

Aigner et al. present in [10] the methodology they developed 
in order to get the hourly load curves of selected domestic 
end-uses for households located in areas next to Los Angeles. 
This analysis technique is in fact the same that Parti & Parti [11] 
developed for the domestic energy consumption analysis: this is 
called Conditional Demand Analysis (often mentioned as CDA). 
Aigner et al. reuse the same mathematical development made in 
[11] (the base CDA) but the authors improved it so that the 
variable which is explained becomes the hourly electric power 
demand. 

As input data, the authors used on the one hand total load 
curves measured at a 15 min time step and on the other hand 
appliance penetration rates for the selected equipments. 
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Fig. 3. Bloc diagram of Aigner et al.’s model [10]. 


CDA model uses 24 regression equations (one equation for 
each hour of a day). Each mathematical relation tries to determine 
the hourly variation of the customers group’s load curve compared 
to the mean daily power demand value.’ The variables which 
explain the power demand are: 


e the desired temperature in the home (two scalar variables): 
they are set with taking into account the outdoor temperature 
of the geographical area’; 

e the size of the dwelling (three scalar variables informing the 
number of rooms); 

è nine dummy variables representing the presence or absence of 
the observed appliances? in a specified household (binary 
variables). 


Running the regression model over the 24h of a day without 
taking into consideration relationships between hours (notably in 
terms of consumption habits) may give incoherent results. That is 
the reason why the authors include in their method some 
simplifications and apply use restrictions. These constitute what 
the authors call a Seemingly Unrelated Regression system. 

Moreover the authors introduce in the mathematical model 
an indicator to differentiate the power demand according to the 
house typology: detached house or flat. The functional architecture 
of Aigner’s model is shown in Fig. 3. 

Aigner et al. present in [10] the results obtained when each 
regression equation is solved separately with three different 
mathematical techniques: 


e ordinary least square method without use restrictions; 
e ordinary least square method with use restrictions; 
e generalized least square method with use restrictions. 


The obtained results with the three calculation ways are in the 
same range. The authors notice that the restriction use is useful in 
order to improve the real aspect of the daily load curves. Yet there 
is no comparison in the paper of the obtained end-uses daily load 
curves with other measurements for each appliance. So the level 
of the regression results are questionable because of this absence 
of validation. 

To conclude with this model, we can say that Aigner et al. 
propose a way to get hourly load curve for some electric 
appliances in a sample of households compared to the direct 


1 This last figure is calculated over an entire month. 

? The three measured samples correspond to so many climatic zones. 

3 Central air conditioner, room air conditioner, water heater, dishwasher, 
washing machine, tumble dryer, range, freezer and swimming pool pump. 


metering of each equipment. This model version can be improved 
if input data concerning the real daily use of the domestic 
appliances is used. In this way, uses of equipments can be defined 
more accurately. Nevertheless the model is limited when con- 
sidering the accuracy of the results. There is no distinction in 
terms of unitary load profile between the same type of electric 
device and the use scenarii between the households. Another 
limit of the model is its disability to supply the load curve of 
appliances which are in a large majority of dwellings: the 
refrigerator is the best example. On the other hand, CDA techni- 
que presents a good ability to identify the load curve of appliances 
with an high energy consumption. 


2.2.2. Bartels et al.'s model 

Bartels et al. present in [12] the load curve model called 
DELMOD they develop on the basis of the CDA technique. In fact, 
CDA is part of its functional procedure but it is not the only one as 
Fig. 4 shows. 

The aim of the authors was to obtain a tool able to simulate the 
impact of various scenarii on the regional power demand of the 
household stock. 

DELMOD is made of two modules: 


e DELMOD B (base module) calculates the load curve of the 
residential sector for a mean working day and a specified 
month. 

e DELMOD W (weather module) quantifies the climatic influence 
on the power demand and modifies the results of the base 
module. 


Results of this model are hourly load curves. 

As input data, this method requires much more information as 
Aigner et al.’s model. In fact, on top of measurements of total load 
curve (400 consumers, 15 min time step) and information con- 
cerning the saturation rates of appliances (which are the same type 
of input data as last model), DELMOD makes use of weather data 
(outdoor temperature knowledge over the load curve measure- 
ment period), socio-economical and technical evolution scenarii 
(expressed through the growth or the decline of economical, 
demographical and social indexes). This feature enables the use 
of DELMOD so as to do demand forecasts. 

In this model, the authors selected nine different appliances.* 
Another variable was included in the model so as to take into 


4 Air conditioner, freezer, automatic defrost fridge, oven, dishwasher, tumble 
dryer, main and/or secondary heater, water heater (off peak and main tariff) and 
pool pump. 
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Fig. 4. Bloc diagram of Bartels et al.’s model [12]. 


consideration the domestic appliances with a large power 
demand. The remaining power corresponds to the in-the-model 
unspecified appliances. 

Some coefficients are used in order to take into account the 
evolution of the technical and/or behavioural aspect which may 
affect the load curve of each end-use. The capacity and the 
efficiency of the appliances are the selected factors. Concerning 
the sociological and economical evolutions three influences are 
used: dwelling size, number of occupants per household and 
income. 

So in order to take into consideration the meteorological 
influences, Bartels et al. introduce in the model the minimal and 
maximal daily outdoor temperature of the considered area (New 
South Wales, Australia). 

From these different sources of information various scenarii 
can be established. All what concerns the non weather influences 
constitutes evolution trends which are constructed in two 
modules of DELMOD B. The weather scenarii are built in DELMOD 
W. There are temperature profiles for the simulation period 
constructed from the measured chronicle on which devia- 
tions around the calculated mean daily temperature are 
implemented. 

Bartels et al. construct the hourly load curves of the selected 
end-uses with help of the CDA technique whose results are 
calibrated using the Model Based Statistical Sampling (MBSS) 
method - developed in [13] - enabling then to obtain better 
results than those from the “simple” Conditional Demand Analy- 
sis. End-uses load curves are then obtained for two types of day: 
with or without evening shopping. 

Bartels et al. simulate the day when the peak load has been 
reached (the “peak day”) and analyze the results of DELMOD. 
Total load curves for customers were calculated after aggregation 
of their different end-uses. Apparently the model enables a good 
representation of the load curve for this day. Other comparisons 
were made for a working day of each month of the year. In this 


case, DELMOD provides not so fine results as for the peak day. 
As explanation for this observation, the authors say that the 
considered end-uses are much more weather dependent in cold 
months. 

DELMOD is a load curve model which can be used for power 
demand forecasts what constitutes a good improvement compar- 
ing with the CDA technique. Various evolution sources can be 
then taken into account. Yet the model suffers the disadvantages 
of the method developed by Aigner et al.: end-uses owned by a 
majority of households cannot be affected with a corresponding 
load curve. Moreover although the hourly power demand for the 
heating system? is calculated by the model, there is no differ- 
entiation in terms of thermal characterization of the envelope of 
the dwellings in the simulated area. Finally, the results of 
DELMOD have been compared on an aggregated basis. That is to 
say that end-use load curves have not been validated individually 
and thus do not have any proven validity for industrial use (e.g. 
impact after retrofiting a dwelling). 


3. Bottom-up models 


3.1. Introduction 

Bottom-up models compute electric demand for a few mod- 
elled households which are either representative of a larger space 
scale or simply target consumers. The unitary results are then 
extrapolated to obtain the electricity consumption for the entire 
studied geographical scale. Here the required input data can be: 


e the individual consumption of the selected domestic appliances; 
e their technical properties; 


5 This constitutes a singularity among the presented load curve models. 


6544 


A. Grandjean et al. / Renewable and Sustainable Energy Reviews 16 (2012) 6539-6565 


of the model 


Daily end-use 
consumption 
per capita 


Random Number 


Equipment 
penetration rate 


Presence at 
home scenarii 


2o 
sx 
2s 
at 
o® 
30 
2g 
éd 
p o 
og 


20 iterations 
of step 1 and 
aggregation 


Building thermal 


model 


Demographic 


data 


Iteration on the 


; households types 


Operational elements 


Generator Technique 


Output data 


End-use load curve for the specified 
human presence scenario 


Iteration on 


the enduses End-uses aggregation 


1 Specific profile 


Load curve for the 
heating system 


Aggregation 


Sum of the Specific daily appliance 


oe PS Iteration on _ profiles and the thermal load curve 
the scenarii Typical household electric 
appliance load profile 
Total domestic load curve 
(1 day) 
Aggregation 
per household 
type 


Regional or national 


load curve 


Fig. 5. Bloc diagram of Yao and Steemers’s model [14]. 


e the geometrical and thermal properties of the modelled 
dwellings; 

weather information; 

electricity bills of households; 

e human behaviour, etc. 


This methodology presents the advantage of a high level of 
details because energy demand can be calculated for a couple of 
households and even for each end-use. With such data, energy 
demand evolutions may be more precisely quantified since the 
participation in the global demand of individual domestic appli- 
ances is known. Energy demand measurements are not absolutely 
necessary and thus it is possible to simulate and to obtain results 
for non-trend scenarii. Among other capabilities, this enables the 
simulation of the impact of technical breaks. Higher precision is 
correlated with higher complexity of the modelling and with the 
requirement of a large amount of input data. 

In this section we present methodologies whose aim is the 
synthesis of residential load curve. With help of accurate data, the 
authors try to build up the load curve of households when 
considering a set of selected appliances and with applying a 
series of assumptions. On top of that, diversity affecting electricity 
consumption has to be generated by the model. Various meth- 
odologies exist to obtain it. To our opinion, some are relatively 
simple, other are well elaborated and complex. In the following, 
the models are classified thanks to the authors’ adopted scale of 
complexity for the generation of diversity. The models appear 
according to a subjective and increasing complexity level of 
diversity. 


3.2. Statistical random model 


Yao and Steemers present in [14] the bottom-up model they 
built up. Called Simple Method of formulating Load Profile 
(SMLP), this method produces end-use load curves with consid- 
eration of various human occupation scenarii. SMLP was estab- 
lished in order to simulate the impact of the integration of 
decentralized energy production, especially with Photovoltaic 
(PV) modules, on the British electricity grid. The aim of this tool 
is notably to predict the peak load curves and to identify the 
network configurations when the system stability can be at fault. 
That is why the authors simulate the worst possible case that is to 
say that all appliances in each dwelling are supposed to run 
through each simulated day. Results of the method are total load 
curves at 1, 5, 15 or 30 min time steps. The model architecture is 
shown in Fig. 5. 

Yao and Steemers make use of different input data: 


e individual DHW profiles; 

per capita daily electricity consumption for each end-use; 

e heating system energy consumption from simulation results of 
Yao et al.’s thermal model [15,16]; 

e statistical results from national population surveys; 

e appliance saturation rates for a series of domestic electric 
devices; 

e constructed occupation scenarii in the dwellings. 


The (15) appliances taken into consideration are listed below: 
TV, VCR, fridge, freezer, fridge-freezer, hob, oven, kettle, micro-wave 
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oven, dishwasher, washing machine, tumble dryer, vacuum cleaner, 
iron and lighting. 

In order to generate the total load curve (all end-uses con- 
sidered) of a selected area and according to the parametrized 
simulation scenario, Yao and Steemers’s model proceeds in three 
steps: 


1. at the household level, generation of one daily load curve for a 
selected appliance and occupants’ presence scenario: the start 
time is selected randomly with help of a Random Number 
Generator Technique (RNGT); 

2. start time diversity is obtained when repeating previous step 
20 times; 

3. aggregation for the specified presence scenario of the daily 
load curves obtained at step 2: the resulting load curve is 
naturally smoother than the individual ones; 

4. repetition of the steps 1-3 for each appliance according to the 
considered presence scenario; 

5. repetition of the steps 1 to 4 for each presence scenario 
considered; 

6. aggregation of previous total daily load curves according to the 
characteristics of the population area (composition of house- 
holds, type of dwellings). 


So as to validate their model, Yao and Steemers used a 
measured load curve coming from the data base of the UK 
Electricity Associate Load Research Team. The chosen load curve 
is a three-person household’s, measured during a winter working 
day. The SMLP simulation to obtain a comparable load curve 
consisted of simulating a 100 household community, to calculate 
its load curve and to divide the result by 100 in order to get the 
load curve of the “mean” household. 
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The used method to give a validity to the model is very 
questionable because of the chosen compared things: individual 
and sample mean load curve are two very different elements. 

To conclude, this model presents the advantage of the simpli- 
city, the availability of the input data and its deeply bottom-up 
character (several aggregation steps are needed). Yet the occupa- 
tion scenarii that are proposed by the authors seems to be very 
simple to represent the reality with a good fit. Moreover the 
introduction of the diversity is done with a pure mathematical 
procedure without considering for instance domestic habits. On 
top of that, the thermal model which is used corresponds to a one 
air-node building modelling which allows a limited daily regula- 
tion of the heating system in a home. There is no differentiation 
and/or no indication in terms of building envelope characteristics 
and appliance efficiency. That is the reason why technical evolu- 
tion scenarii cannot be simulated with this model. 


3.3. Probabilistic empirical models 


3.3.1. Stokes’s model 

Stokes developed a model applicable to the domestic sector 
and presented in [17]. It is able to generate three types of 
aggregated load curves 


1. 30 min demand profile for an average household (1st level), 
2. 30 min demand profile for a specified dwelling (2nd level) and 
3. 1 min demand profile for a selected domestic customer (3rd level) 


with a fine spatial resolution that is the individual domestic 
appliance. 

To get this gradual precision, the model is structured into three 
main modules; results of the level n serve as inputs of the level 
n+1. Fig. 6 schematically illustrates this architecture. 


Total load curves 
1200 clients, 30 min time step 
with individual annual peak 
load estimation 


End-use load curves 
175 households, 9 end-uses 


Operational elements 


of the model 


Data statistical analysis 


Sinusoidal end-uses 
modelling 
(active power) 


Output data 


Average household total 
30 min load curve 


magnitude s 


Heating and DHW data 


Socioeconomic data 


Appliances data 


- nominal power 


random 
parameter r 
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- unitary load cycle 


std-deviation o 


Reactive power calculation 
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Krightiow Kitestyle 


Probabilistic triggering p of 
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Fig. 6. Bloc diagram of Stokes’s model [17]. 
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Because of the lack of 1 min measured load curves which can 
be used as (calibration) input data and in order to avoid modelling 
precision discrepancies between appliances, Stokes built up a 
unique, clear and homogeneous modelling framework that is 
applied to all end-uses. This way, the level of details of the 
modelling is the same whatever the domestic device. On top of 
that, Stokes’s model calculates the load factor of each equipment 
at each time step enabling then to take into account the reactive 
power demand. According to the author, this part of the power is 
often neglected in load curve models despite its interest as for 
instance for load flow analysis purposes. 

In the following, the input data for each level, each computing 
phase with focus on the assumed simplifications, the results of 
each level and finally the validation of them will be considered. 

Input data: In order to calculate the 30 min demand profile for 
an average household, the 1st level of Stokes’s model uses 30 min 
load curves that have been measured in a sample of 1200 English 
households. Moreover data at the device level (individual load 
curves for nine appliances) and knowledge about production and 
stocking of DHW: this set is differentiated in terms of day-type 
(three modalities) and tariff (two options). On top of that, 
information about socio-economic characteristics of the house- 
holds on the one hand and concerning the dwelling on the other, 
was available through answers of questionnaires. With regard to 
the national average, different discrepancies have been identified 
in the sample by the author. 

On top of 30 min demand profile for an average household, the 
2nd level of the model, whose aim is the calculation of 30 min 
demand profile for a specified dwelling, receives information 
about the penetration rate for each appliance, the occupation of 
the households (number of occupants according to the dwelling 
area and presence at home scenarii) and an estimation of its 
annual maximal power demand value. The type of dwelling, its 
surface area, its orientation, etc. are other diversity elements that 
have been considered out of Stokes’s work because of the lack of 
so precise data. 

So as to step down at the 1 min time resolution and to provide 
corresponding load curves for a specified household (load curve 
that is derived from the 30min load curve calculated by the 
2nd level), the 3rd level gets information concerning the habits 
of appliance use (e.g. duration and frequency of the cycles), 
measured unitary load cycles and estimations of maximal power 
demand per appliance. According to these data, four parameters 
are defined to characterize the unitary load cycles at the 1 min 
time resolution: scale of demand, event duration, use frequency 
and timing. 

Modelling details: Stokes’s purpose is the modelling of the 
domestic load curve for average British dwellings. Her model has 
to be able to provide demand profiles for contemporary situation 
as future one: so load curves must be adjustable. That is the 
reason why she decided to assume following statement: 


Variations that affect end-use peak demand and its annual 
consumption are identical in the time. 


Concretely, reduced (normalized) power demand, P,.;, for each 
end-use has been calculated thanks to the Load Research Group® 
(LRG) data set according to Eq. (1) 


P 2 
Pmodel,i => = Py i a) 


P MAXq,j 


where i is a specified end-use, P noaei is the end-use power data 
used in the model, Pmes,i is the measured 30 min power demand 
for end-use i and Pmax,; is its annual maximal power demand. 


6 See [17] p. 6. 


From the LRG data set, Stokes tried to build up a mathematical 
approximation for 30 min power demand fluctuations on the 
year. According to her research, the best compromise between 
accuracy and simplicity is the sinusoidal modelling of the 30 min 
normalized power demand — Pyinusoidai — that is reported in Eq. (2). 


. N 
Pinusoidal =s. sın (27 : x) -¢] +k+r 2) 
y 


s represents the sinus function amplitude, @ its phase, Ng 
corresponds to the considered day’s number in a year with N, 
days. k is a constant and r is a random number generated with an 
end-use specific Laplace-Gauss density function [0,¢sinusoiaatl- 

All parameters derive from the LRG data set. They intrinsically 
include weather influences: that is the reason why Stokes’s model 
does not make use of “traditional” weather file as input data. 

Moreover because of the independence between the Laplace- 
Gauss density functions, linked end-uses are not considered in 
Stokes’s model. 

For certain end-uses (e.g. lighting), minimal and maximal 
power demand levels have been defined and are applied when 
Eq. (2) is not relevant or gives negative values. In these cases, 
30min reduced power demand P, is calculated according to 
Eq. (3). 


l Prinin +T min if P, < Prin 
r= 


if P, > Pr (3) 


Pr nax + Vmax 


Tmin (respectively Tmax) is a random number generated with a 
[0,0 min] (respectively [0,¢max]) Laplace-Gauss distribution. 

Long term consumption trends have been provided by the 
Market transformation Program [18] for each end-use i. From 
these, trend estimations of end-use maximal power demand 
Prrend.max,i Nave been calculated enabling then the scaling up of 
the normalized demand (P,,;) generated by Eqs. (2) and (3). Thus 
households’ group average 30 min (non reduced) power demand 
— Pgroup,i — is calculated with help of relation (4). 


Peroup,i = Pii ý Prrend,max,i (4) 

In her dissertation [17], Stokes adapted the previous modelling 
approach for various individual end-uses listed in Table 1. Further 
details about their modelling and the taken assumptions can be 
directly read in the referenced document. 

The results of this level have been validated on measures 
conducted by the LRG. Comparisons in terms of 30 min aggre- 
gated consumption (all end-uses considered) showed the good 
ability of the model in terms of consumed energy. 


Table 1 
End-use families and appliances considered in Stokes’s model [17]. 


End-use family End-uses considered 


Domestic cold Fridge 
Freezer 


Combined fridge/freezer 


Heating Storage heater 
Domestic hot water Immersion heater 
Cooking Hot plate 

Oven 

Kettle 


Micro-wave oven 

Washing-machine 

Tumble dryer 

Combined washing-machine/tumble-dryer 
Dish-washer 

“Standard” bulb 

Other appliances 


Wet appliances 


Lighting 
Miscellaneous 
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The second level of Stokes’s model introduces some influences 
responsible for diversity in order to produce the 30 min power 
demand load curve for a specified household. 

Income - in relation to the English mean salary - and lifestyle 
are two elements that influence the use of electricity because 
household’s type of dwelling and corresponding room number are 
correlated to them: that is why they were integrated into the 
model. However they are only applied for lighting and miscella- 
neous appliances on the one hand, they influence the equipment 
set of the households on the other. 

According to the dwelling surface area, that is linked to the 
number of rooms, the household size (number of members) is 
calculated. Then household’s consumption for each end-use is 
adjusted according to its size notably with the help of results from 
Boardman et al. [19] study. 

However the presence at home daily profiles and household’s 
occupants’ age are two variation sources that are not taken into 
account in Stokes’s model. Nevertheless correlations between 
lifestyles, penetration rates and daily appliance use were identi- 
fied in [20] and applied when possible. 

In order to calculate 30 min (Half Hourly) load curve for a kind 
of household (not a precise one) and an end-use i, Py spei, the 
second level of Stokes’s model declines Eq. (5) 


Punt spe,i = Putt group,i . Kownership,i i Koccupation,i (5) 


where Pu,group,i is the calculated half hourly power demand for 
the end-use i into a group of customers, Kownersnip,i is a boolean 
factor indicating the presence (or absence) of the end-use i in the 
considered household and Koccupationi is a Scalar factor that modi- 
fies the value of the calculated power demand according to the 
size of the regarded household. 

For some end-uses i, another coefficient (kjgestyie,) multiplies 
previous equation. It is aimed to take into account the lifestyle 
(more or less the social status) of the simulated household. For 
the domestic hot water, Kiifestyle is replaced by a coefficient (four 
modes: very low, low, average and high) whose mode is randomly 
selected that adjusts the water usage when simulating a sample 
of households. 

The output of the second level of Stokes’s model is calculated 
when adding all the previously determined 30 min load curves for 
each end-use. However this level does not provide results for a 
specified household because all dwellings with the same char- 
acteristics (household size, lifestyle, domestic equipment set, etc.) 
will be affected with the same load curves. What introduces 
variations between them is the third level of the model. It takes 
the previous calculated power demand profiles as input and 
considers them as probability distributions to start the end-uses 
within a day. 

In the third level of her model, Stokes tries to explain in detail 
the diversity of the domestic energy consumption in order to 
provide the calculation of the 1 min domestic load curve for a 
specified household. According to Stokes, power demand cannot 
be modelled without consideration of the chance influence. That 
is the reason why the end-use launching is determined by 
comparison between a generated [0; 1] random number and a 
constructed ratio p calculated as follows: 


Pur spe,i(hh,d) 


p(event occurring) = 
P HH,app,i 


(6) 
where p(event occurring) is the probability that an event (the 
launching of an end-use) occurs, Pypspei(hh,day) is the half hourly 
power demand value for the end-use i calculated for a specific 
household at the day d and the half-hour hh and Py app,i is the 


power demand required for a “standard” appliance app that 
serves the end-use i. 


The more the value of the power demand calculated at the 
second level is high, the more the probability p to start at least 
one appliance is important. On top of that, the model is able 
to start several equipments that are related to the same end-use 
(e.g. lighting) on the one hand, it takes into account the different 
power demand levels when relevant on the other hand (e.g. hot 
plates). 

Concerning the time start selection of the end-uses, they 
are chosen with considering individually each half hour. Apart 
from equipments whose functioning cycle(s) is (are) longer 
than 30 min, the appliances start times are chosen so that the 
end-use required energy demand is satisfied into the current 
half-hour. 

The end-use starting frequency and the duration of each cycle 
are two elements that make the model adaptable. Other assump- 
tions have been made by Stokes for each kind of end-use. 

To validate her model, the author makes use of two data sets: 


e measures on a bungalow district at the 5 min resolution; 
e measures on a sample of 13 households at the 1 min time step 


and considers different comparison criteria: energy demand, 
maximal and mean power demand value (power factor), daily 
profiles, load distribution density function. 

On the first sample, Stokes’s model succeeds in estimating the 
energy demand on the simulation period. However, it seems to be 
the result of various overestimations and underestimates depend- 
ing on the considered end-use. In terms of power demand, 
generated daily profiles are relevant but measured maximal 
power demand value is always greater than the simulated one. 
According to the author, it seems to be due to the too low 
generated diversity. 

On the second sample, the same conclusions can be said but on 
top of that the model gives a general underestimate of the energy 
consumption and is quite unable to predict the morning and 
evening peak loads of the domestic customers. 

To conclude, Stokes established a very detailed model split 
into three levels. She adopted an homogeneous modelling frame- 
work that makes easier its reuse and modification. Differentiation 
between end-uses and appliances is carefully taken into consid- 
eration. Moreover, the author tried to implement in her model a 
large set of power demand determinants. 

However as we said in the validation paragraph, the genera- 
tion of diversity reaches its limits because of the use of random 
coefficients that do not really represent the existing end-use 
relations. On top of that, the 1 min load curve result for a specified 
household is partly dependent on the group of measures used for 
the calculation (because of the load curve “break down” process). 
Moreover Stokes’s model does not explicitly take into account the 
daily occupation scenarii of the occupants for the end-uses 
starting algorithm (it is supposed to be integrated in the input 
data). Finally, the heating power demand and associated energy 
consumption is only dependent on the dwelling size area: yet the 
building characteristics, especially the insulation level of the 
house, plays a determinant role for the thermal appliances. 


3.3.2. Paatero and Lund’s model 

In [21], Paatero and Lund explain the model they have set up. 
Its aim is to provide 1 h load curve for domestic customers (from 
few dwellings to a relative large set—1000 clients). Its architec- 
ture is represented in Fig. 7. 

The authors made use of two main data set as inputs. In fact, 
there are 1h power demand measurements that were led on 
Finnish domestic customers living in flats (respectively during 
1 year on 702 households and during 143 days on 1082 
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Fig. 7. Bloc diagram of Paatero and Lund’s model [21]. 


dwellings). Because of the absence in the measured apartments of 
electric heaters, coolers and immersion heaters for domestic hot 
water supply, these end-uses are not simulated by the model. On 
top of that, Paatero and Lund used various studies to construct a 
data base of unitary load cycles and standby power consumption 
for the simulated appliances. Moreover, they gather national 
statistics in order to build up the domestic equipment sets. 

In order to develop the model, Paatero and Lund studied a 
series of load curve measures. This way they noticed that power 
demand fluctuates according to different periods especially 
the season, the day and the hour. Exogenous influences like 
outdoor temperature and daylight availability can be theoretically 
responsible for them. However in this dataset, these influences 
cannot explain such fluctuations because of the absence of 
thermal appliances. Domestic habits and social behaviour are 
responsible for this variability but they did not have been taken 
into account because it constitutes a too high level of details 
according to the authors. Moreover Paatero and Lund identify as 
well as Capasso et al. the deeply interactions between the 
presence at home of the households’ members and conducting 
domestic activities. 

Paatero and Lund constructed a two-level model: the first level 
defines both the domestic equipments set for each household and 
the general load fluctuation trend. The second level simulates 
each appliance in each household with the help of end-use 
starting probabilities. 

To generate diversity, the first level calculates the density 
function of the domestic daily consumption with help of sample 
data. The authors have reduced the previous distribution in order 
to get a dimensionless curve that is centered on the daily mean 
consumption. Then for each simulated day, a value is selected 
according to the distribution: it is called Psociaı (social random 
factor) because it is supposed to capture the social diversity of the 
demand and it is applied identically for each simulated house- 
hold. The equipment set in each dwelling is constructed according 
to the ownership statistics. 


The second level of the model provides the hourly load curve 
for each equipment with considering each household separately. 
To compute this, two elements are required: the unitary load 
cycle of each appliance and its starting probability Pstart calculated 
according to Eq. (7). 


Pstart (1, we,At,Gap,hr,d) = Pseason(l, We) a Prr(i,hr,d) - fd) 
. Pstep(At) : Psocial(F ap) (7) 


with: 


i is an equipment or a group of equipments, 

hr is the hour of the day, 

d is the day of the week, 

we is the week of the year, 

Psociat is the social random factor whose standard deviation is 

Cap (ap is for apartment), 

© Dseason is the seasonal factor (models the weather influences), 

© Pn, is the hourly factor (that takes into account the changes in 
terms of activity levels during the day), 

e fis the mean daily starting frequency for each end-use i, 

e At is the working time step and 

© Pstep is a scaling factor that adjusts the probability values 

according to At. 


According to the sign of the comparison between Pstart and a 
[0; 1] generated random number, the appliance i is started or not. 
If the equipment is effectively launched, its load cycle is added to 
the currently calculated household total load curve. Apart from 
“cumulative” appliances (such as light bulbs) the calculation of 
Pstart is stopped for the equipments during their functioning. As 
soon as the end of the cycle is reached, the appliance returns then 
in off-mode and the starting probability can be evaluated again. 

Concerning the standby consumption, it is taken into consid- 
eration with adding a constant load at each time step. 

As a exercise, Paatero and Lund tested their model to calculate 
the yearly energy consumption of a 10000 households sample 
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living in flats. It is not a validation of the load curve reconstitution 
method because no comparison is carried out on the generated 
power demand profiles. They only conclude that their model 
underestimates the mean daily energy consumption and that the 
distribution of the calculated daily energy use is more spread as 
the measured data. 

Paatero and Lund developed a bottom-up model that is able to 
provide load curve that are used for instance to quantify the 
impact of DSM measures. The strength of this method is that it 
relies on public data and its computational process is relatively 
simple. 

Yet the model is limited in different ways. First of all, because 
of the used input data, the method is only able to calculate load 
curve for households living in flats: the data sample is highly non 
representative of the entire housing stock. Then the model cannot 
really step back from the past reality because the random 
influences are constructed on the input data set. On top of that, 
the model is not capable of evaluating the power demand impact 
that follows the introduction of new end-uses, appliances or 
behaviour. 


3.4. Time of use based models 


Time of use data - whose concept was originally developed by 
Pratt et al. in [22] - generally correspond to daily frequency 
curves of use for each end-use/domestic activity. In order to 
construct this kind of information two main ways are available: 


e a monitoring campaign is realized in a sample of households 
and during the study time each household’s member has to 
write down in diaries the conducted domestic activities at the 
time step resolution that is chosen for the campaign, 

e the other way is to focus the study on different general 
domestic activities (breakfast, lunch and dinner time, get up 
and bedtime time, daily transport duration, etc.) and to 
reconstitute an incomplete diary with mean times of each 
activity with a density function around these moments. 


Then the compilation of the results enables one to obtain 
frequency curves for the studied activities by comparing for each 
time step the number of households where a specified activity is 
being conducted with the total number of studied households. 
The segmentation of the results can be done for each end-use 
with the same methodology if the study has been conducted on 
the domestic equipment basis. Moreover it is possible to get 
curves for different types of day (it is often divided between 
weekday and weekend day). This concept enables the modelling 
of domestic activities and it provides realistic sequences of 
electric equipments use in a sample of households. 


3.4.1. Walker and Pokoski’s model 

Walker and Pokoski [23] are the first authors which use the 
concept of time of use with the aim of constructing a load curve 
model. More precisely they are the first who take into account 
human behaviour to reconstruct load curves for a specified 
household or a set of different dwellings. The model was devel- 
oped so as to help the planning of new power plants and it gives 
total load curves at 15 min time step. 

The basis of the human behaviour modelling is the use of two 
kinds of probability functions respectively called availability and 
proclivity which will be then generalized in the model of Capasso 
et al. (cf. Section 3.4.2). The first one gives the household’s 
occupants probability to be at home and so “available” to use 
an electric domestic appliance. Concretely it is implemented in 
the model with 96 x 1 daily vectors (the time resolution is 15 min 


- 96 x 15=1440 min=24 h) containing 0 and 1 indicating respec- 
tively on the one hand the absence or the asleep state of a person 
and, on the other hand, the presence at home of a specified 
household member in an awake state. Theoretically in order to 
represent the exact reality, the input data required for this 
function is one occupation diary per studied person. Consequently 
these data-gathering can be very expensive. That is the reason 
why Walker and Pokoski only make use of one presence scenario 
and they improve it with variations sources thanks to different 
density functions implemented at the start times of the main 
activities. In this way they get various “mathematically con- 
structed” occupation daily scenarii. Yet the authors had to take 
into account the different ways people live: typically the time to 
go to work and to come back from work generally in the evening. 
Delay probabilities were defined in accordance with the esti- 
mated proportion of each “way of life” in the simulated area. 
All things considered, the availability function is constructed 
with help of 13 sub-functions, each one gives a unitary kind of 
information. Moreover Walker and Pokoski take into considera- 
tion the absence from home because of non professional reasons 
(hobbies, shopping, etc.): they model this influence with the help 
of Paway, a probability function which makes dependent on it the 
availability function. Paway is compared with a random number 
generated with a [0; 1] uniform density function. Two ways are 
then adopted in order to calculate the number of occupants at 
home for each time step according to whether Paway is bigger or 
smaller than the generated random number. 

The availability function gives not enough information in order 
to reconstruct the load curve of a set of households. That is the 
reason why the authors make use of various proclivity functions 
which model domestic habits and more precisely time periods 
when they are likely to be conducted.’ The proclivity for a certain 
activity is equivalent to the tendency to do it during a day. These 
preferential time periods come from standards of living, conven- 
tional ways of life or even people habits. In fact there are two 
kinds of proclivity functions: 


1. the first one indicates the likelihood for the use of a specific 
individual electric appliance; 

2. the second type considers a domestic activity which may 
require the use of various equipments. 


In order to construct the first proclivity function (the so called 
basic proclivity function), Walker and Pokoski make use of 
different results coming from various surveys whose aim was 
the construction of a daily diary for a selection of domestic 
appliances so as to know the periods when they are likely to be 
used and the weekly frequency of turn on events. 

The second type of proclivity functions represents more gen- 
eral activities such as the daily mealtimes. Potentially this activity 
involves the use of a few domestic equipments such as the fridge, 
the range and/or the water heater because of the very possible 
draw of water for this activity. With help of results from various 
studies, Walker and Pokoski define density probability functions 
for each mealtime of the day. This way each period in the day is 
affected with a probability that meal can be taken during this 
time period. Such probabilities can be defined for other activities. 

With the help of the availability and the proclivity functions, 
the authors have the knowledge concerning the moments when 
the electric appliances are turned on during the simulation time. 
To get the corresponding load curve for each end-use they have to 
“stick” on the identified on/off events, unitary load curve patterns 
that are specific for each appliance. In the case of the appliance 


7 Without excluding their respective probability at any other time in a day. 
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Fig. 8. Bloc diagram of Walker and Pokoski’s model [23]. 


which run without the need of a human intervention (for instance, 
they are turned on or off according to the information which is 
provided by sensors), Walker and Pokoski use what they call 
“normal-cycle” functions. It is notably the case of the fridge which 
runs all the time with a cycling on/off alternation. 

The architecture of the model is presented in Fig. 8. 

Fig. 8 shows the chosen aggregation method. In fact, the 
simulation tool provides individual end-use load curves (one curve 
per electric appliance). These raw results are then aggregated 
according to two steps: 


1. first a summation is done at the household level achieving 
then the total load curve for a specified dwelling; 

2. secondly the total load curves are aggregated in order to get 
the global load curve of the simulated area. 


This modular architecture enables one to consider the influence 
of the power demand of an individual appliance with regard to the 
total demand of the household. Another advantage of this technique 
is the possibility to add a new appliance in the model without 
reconstructing its whole structure and with conserving the calcula- 
tion method. Moreover the simulation provides the load curve of an 
area which is not calculated as an individual household’s mean load 
curve multiplied by the number of simulated dwellings. Thus real 
behaviour with potentially huge differences between them can be 
simulated. Furthermore it gives the opportunity to simulate a lot of 
dwellings with different sets of domestic equipments. 

In [23] only the water heater is presented in detail but the 
entire model is described in [24] and at the date of the publication 
the model did not contain thermal load appliance (such as 
heater). However Walker and Pokoski did not consider this 
feature much more complicated than the water heater because, 
according to the authors, it only requires another function which 
takes into account, on top of the availability and proclivity 


information, the in- and out- heat flows affecting the dwelling 
and economic considerations concerning the occupants. 

To validate their model the authors carried out two simula- 
tions. The first one had the aim to obtain the daily load curve for 
two single-family homes. The second test concerned two groups 
of residences. Simulation results were compared with data from a 
load curve measurement campaign. The chosen criteria for the 
comparison was the Normalized Variation Factor (NVF) defined as 
follows: 


ie 1 sn 

data 1Pese(t)—Pmes()” z Nea 1 Pese(t)—Pmes(t))” 
pH 2 = 1 ‘ 2 
Pes G . St -1 Presl) 


with Pes{t) the estimated power demand at time t, Pmes(t) the 
measured power demand at the same time and n the number of 
simulation time steps. Conceptually the NVF corresponds to a 
“pseudo-variance” of the estimated power demand calculated 
with reference to the measured load curve that is normalized by 
the square of the mean measured power demand value. So NVF 
represents the dispersion of the estimated power values com- 
pared to the measured ones. NVF is reused by Capasso et al. [25] 
and Widén et al. [26]. 

With this indicator, the authors tried to know if their results 
are included in the fluctuation range existing between two daily 
measured load curves. In fact, it is possible to calculate the NVF for 
the same household (or group of households) between different 
days or for the same day and between different households. They 
got good NVF values (i. e. small figures) confirming the fact that 
the simulated load curves fit the measured ones. 


NVF = 


(8) 


3.4.2. Capasso et al.’s model 
Capasso et al. present in [25] the model they developed which 
is called ARGOS. Its aim is to reconstruct the load curve for a 
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specified household up to a set of dwellings, in a defined area, in 
order to quantify the impacts of DSM measures. These are taken 
namely so as to provide a solution for one of the major problems 
affecting the Italian electricity network: the winter peak load 
occurring during worked days. Another purpose of this metho- 
dology is to be a tool for optimizing the production of the power 
plants or even to “drive” the planning of new energy supply sites. 
The main characteristic of this approach is the fact that it takes 
into account the major interactions between the electricity 
energy demand and the behavioural factors. The results are load 
curves with a 15 min time resolution. 

From results of various surveys, the authors constructed some 
probability functions. They notably make use of the Walker and 
Pokoski’s availability and proclivity functions enabling then to 
obtain mean profiles for the basic proclivity function. 

ARGOS simulates the load curve of a residential customer. This 
is built up with a set of domestic appliances and a specified 
number of occupants. According to Capasso et al. the real issue 
for calculating domestic load curves is the modelling of the inter- 
actions between the set of equipments and the household’s 
members. 

So as to set a realistic modelling frame, rules and constraints 
are programmed in ARGOS under a probabilistic form: there are 
behavioural and engineering functions. 

On the one hand, the first ones define the appliances set for 
each household according to its characteristics. On the other 
hand, they establish the specified household corresponding time 
structure for the use of the equipments. That is to say that the 
domestic activities, which are grouped together in four main 
families (cooking, housework, leisure and personal hygiene), are 
not randomly conducted by the occupants through a day, the 
starting of the equipments too. On the contrary, they are both 
spread in the simulated day thanks to time of use density 


functions. These are calibrated so that simulations lead to 
appliances daily consumption that are included within a con- 
fidence interval around values coming from measurements. 
Simultaneously each household member is characterized with a 
presence at home scenario and a specific proclivity for each 
domestic activity. 

The engineering functions give various information: first they 
precisely characterize the functioning of each domestic equip- 
ment according to its properties, secondly they set each house- 
hold subscribed maximal electric power and finally they set the 
technological penetration in the simulated population. 

The previous described data and probability functions (beha- 
vioural and engineering ones) are stored in a file which actually 
contains all required inputs of the model. 

From the parameters that describe the scenario, a Monte Carlo 
extraction is carried out to choose the simulated type of house- 
hold. The corresponding total load curve is calculated in ARGOS 
by addition of each end-use specific load curve. 

The total load curve for a complete area is simply calculated by 
addition of the individual consumer’s load curves. The architec- 
ture of ARGOS is shown in Fig. 9. 

As model validation task, Capasso et al. tried to reconstruct the 
aggregated load curve of 180 apartments which are representa- 
tive of the Italian dwellings located in the main cities’ suburbs. 
Among them, only 95 households have been kept as comparison 
references because they answered a questionnaire whose results 
built up the ARGOS input data file. 

As we previously said, this model was established especially to 
predict the residential load curve during peak load days. That is 
the reason why, the authors focused the comparison results on 
restricted time periods. For these critical hours, modelled and 
measured load curves are very similar. For other infra-day 
periods, comparisons are less favourable. 


- Recognition label 
- Operation mode 
- Power demand limits 


- Duty cycle duration 

- Daily usage profile and probability 
- Unitary energy consumption 

- Penetration rate 

- Required human ressources 

- Deferrability 

- Unitary load cycles 


Appliances (specific elec.) 


- Occupant type 

- Income 

- Daily presence profile 
- Proclivity 


Occupant 


- Occupants types 
and number 


- Mean electricity 
consumption 

- Electric appliances 
ownership 


SJaquiaW S,pjoyasnoH 


Household 


selection 


Other input parameters : 
simulated day, geography, weather... 


Operational elements 


of the model 


Modelling of the er 
interactions 


Behavioural and 
technological functions 


- Number of identical Variable 
households client i 


Output data 


Regional or national 
load curve 


Aggregation of n 
clients In the area 


! n simulations 


Ww 
co 
23 Domestic client total 
Os | 
D ge oad curve 
Monte Carlo 5 S (except from heating and DHW) 
extraction =1s 


ee ee te 


End-uses load 
curves aggregation 


End-uses load curves 
(except from heating and DHW) 


pne P eee en een een e one eneeneco eee nocoecoscoscoscoscos 


Fig. 9. Bloc diagram of Capasso et al.’s model [25]. 
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In a second step, the comparison period was extended to a 
complete simulated day. For this task, Capasso et al. estimated and 
compared the NVF between measured and modelled load curves on 
the one hand and between two consecutive daily measured load 
curves on the other. Results satisfied the confidence interval 
criterion indicating the good approximation ability of the model. 

ARGOS constitutes a noticeable improvement of the Walker 
and Pokoski’s model and in the same time a comparison reference 
for the models developed after it, especially for the human 
behaviour modelling. Finally, this is notably well reproduced 
because of the use of real elements instead of standardized 
occupation profiles as in Yao and Steemers’s model. Another 
advantage of ARGOS is its ability to provide the load curve results 
under various aggregation levels enabling then the calculation at 
different scales of load curve specific coefficients. Moreover, new 
appliances can be implemented in the model, rules and con- 
straints are changed and adapted by the user. 

Yet, these possibilities require a lot of input data that are 
difficult to get at a large geographic scale (i. e. the simulation 
of the national domestic load curve seems to be illusory). On 
top of that, the needed information for the elaboration of the 
behavioural functions have been systematically adapted into 
Laplace-Gauss distribution. This can lead to non optimal adjust- 
ments in all cases. Some significant distinctions in setting the 
parameters of the simulation have not been taken into account: 
worked day/not worked day, winter/summer, etc. Nothing about 
building diversity (for instance thermal characteristics) has been 
implemented in ARGOS because of the relative insignificance of 
electric heating in Italy. Finally, the influence of the electricity 
tariffs is not supported by this model. 

After the first model implementation in FORTRAN, research 
has been resumed and continued by Prudenzi, Falvo and Silvestri. 
This work leads to a new software called DELOS (see [27,28]). 
However no significant evolutions have been added in this model 
compared with the initial one (ARGOS): in particular the concepts 
remain the same. 


3.4.3. Armstrong et al.’s model 

In the framework of Annex 42 of IEA (International Energy 
Agency) “The Simulation of Building-Integrated Fuel Cell and Other 
Cogeneration Systems (COGEN-SIM)”, Armstrong et al. [29] estab- 
lished a model which reconstructs the load curve of domestic 
specific electricity appliances with a 5 min time resolution. The 
aim of the modelling task is to obtain various typical consumption 
profiles for the selected appliances so as to determine the way micro 
combined heat and power (u-CHP) can be integrated in Canadian 
dwellings and how they can provide the required power. 

In this frame, three target households have been simulated. 
They differ in terms of dwelling surface area and intensity of 
energy consumption. 

Armstrong et al. make use of information and notably: 


e estimated target household annual energy consumption for 
each selected appliance® thanks to three sets of use factors® 
that compare actual with national average appliance energy 
consumption, 

e domestic equipment set with their characteristics (inclusive 
penetration rates!°), 


8 Fridge, freezer, dish-washer, washing machine, tumble dryer, range, lighting 
and other appliances. 

9 Non validated coefficients. 

10 Because of the choice of the target dwellings, the authors adapted the national 
penetration rates. Concretely they systematically assigned for each household some 
selected appliances (fridge, freezer, dish-washer, washing machine, range and 
tumble-dryer) and they applied target household specific hypothesis. 


e and occupants’ use profiles. 


This wide set of information was established thanks to 
surveys, statistical and appliances manufacturers data; it is 
representative of the situation of the Canadian domestic appli- 
ances stock in 2003. 

On top of that, information about unitary load cycles, cycles 
duration, estimated month equipment use and appliances typical 
nominal power have been used by the authors. 

Moreover time of use data have been implemented in the 
model. Because of the lack of these kind of information, Arm- 
strong et al. reused the TOU profiles defined by Pratt et al. [22] for 
some appliances (range, dish-washer, washing machine) without 
any modification. For other devices, they performed some adap- 
tations on the available profiles. 

The load curve construction methodology, which is schemati- 
cally illustrated in Fig. 10, relies on simple principles and 
assumptions regarding on the one hand the timing and on the 
other hand the magnitude of the power demand. 

The management of the appliances starting is performed with 
the help of the TOU curves and an appliance specific chance factor 
c that provides the estimation of the hourly probability p of a start 
event according to the equation: 


i (9) 


where f is the hourly fraction of the daily use (hourly TOU value). 
For each household’s appliance, c value was set by iteration so as 
to meet the target consumption via the annual cycles number. 

With regard to the magnitude of the load curve construction 
and despite the individual simulation of each device, less model- 
ling methods have been developed than the number of different 
appliances. So it means that a same method was applied to 
several devices whose functioning showed similarities. 

In fact, simplified modelling approaches have been chosen by 
the authors: 


e assignment of a constant power demand during the functioning 
of washing and cooking appliances; 

e the power demand of domestic cold appliances is directly 
constructed from measures that are adapted by the authors to 
match the consumption targets; 

e the load curve of lighting and other appliances is constructed 
from the TOU curves, a duration use that is randomly selected 
between chosen lower and upper bounds and nominal power 
values; 

e a 65W base load represents the stand-by power demand of 
the domestic equipments in the household. 


Then Armstrong et al. proceed with an aggregation step to get 
the daily load curve for a specified household: the eight previous 
calculated loads (one per appliance) are simply added up. The 
yearly load curve for each dwelling type is obtained when 
performing multiple simulations. 

In order to evaluate the goodness of fit of their model, the 
authors carried out comparisons between simulation results and 
load curve measurements. These have been collected by Hydro 
Quebec in the 1990’s during a campaign focused on electric 
heated detached houses. With this campaign total, heating and 
DWH loads have been simultaneously measured at a 15 min time 
resolution. By difference, specific electricity load curves were 
obtained. Four houses from this campaign have been selected to 
proceed the evaluation task. 

First, the visual comparison showed the good ability of the 
model to be near the reality in terms of peak loads, mean loads 
and annual energy consumptions. 
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Fig. 10. Bloc diagram of Armstrong et al.’s model [29]. 


Secondly, Armstrong et al. notice that the model produces less 
repetitive results than what is seen in the measurements. For the 
authors it does not constitute a weakness of their model since the 
generated loads will be used as inputs for the simulation of 
L-CHP: variable demand profiles lead to stress on the CHP that are 
more difficult to comply with. 

Thirdly a statistical comparison (study of the repartition of the 
values in 100 W ranges) has been proceeded leading to two main 
observations: 


e the selected power level for the modelling of stand-by modes 
(65 W) is not optimal; 

e the second observation is more trivial: in measures, two 
100 W-ranges of power demand are empty what is not the 
case of simulated results. This can be explained by the 
domestic appliances sets in the selected households. 


This load curve reconstitution tool shows a series of advan- 
tages: first its relative simplicity, secondly the reasonable 
required amount of input data compared with the other models 
and thirdly its ability to generate various specific electricity 
demand profiles. Yet some weaknesses can be found: the used 
TOU curves probably have few validity because they are not 
recent (1989): domestic habits can evolve significantly. Moreover, 
apart from the used load curves for refrigerators and dishwashers 
that come from measures, the other unitary load cycles are built 
up when assuming a constant power demand during the func- 
tioning what is a real simplification of the reality. Finally the 
model does not make the difference between the types of the 
days (weekday, weekend days for instance) and the influence of 
the weather only concerns the lighting and it is very simplified in 
this case. 


3.4.4, Widén et al.’s model 

Widén et al. expose in [26] the load curve generation model 
they built up. This TOU data based model makes use of simplified 
end-use to power demand conversion schemes, natural daylight 
density function and water consumption; moreover it defines and 
simulates “mean”'' domestic equipments. 

The main purpose of this model is to generate electric load 
curves for domestic appliances that will then serve as inputs of 
other simulation tools. Among them are modelling programs of 
small scale decentralized electricity generation systems or 
domestic water production equipments based on solar energy. 
Another purpose of this simulation tool consists of studying the 
modifications on the domestic load curve in a prospective point of 
view when considering the changes in terms of behaviour and 
energy efficiency of the appliances. 

In the authors’ opinion, a load curve generation model can be 
classified according to two criteria: its spatial and temporal 
resolution. 

The spatial resolution corresponds to the smallest simulation 
level where a tool is able to provide a curve of need (electricity, 
heating, water, etc.): from a unitary domestic equipment to a 
household type to the national residential field. The temporal 
resolution of a model is linked with the time step results are 
returned: from seconds to hour. 

In terms of input data, the main requirements of this model 
are time of use information with a 5 min time resolution that 
have been collected during a series of campaigns. Concretely, 
it is constituted by empirical presence at home scenarii and 
domestic activity sequences for a weekday and a weekend day. 
During them, each occupant of the selected households had to 


11 Tn terms of performance. 
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write down his activities and the corresponding start and stop 
times of them. That notably includes the instants when getting 
up and to go to the bed and when moving out and to home. 
Thus the model provides a load curve for each household’s 
member and can construct load curves for a mean domestic 
customer. 

On top of time of use information, end-use load curves 
measurements and domestic water profiles have been used for 
the elaboration of this model. 

In the whole, five surveys have been used as input data sources, 
comparison elements data bases and validation references. 

The original characteristic of Widén et al.’s model, whose 
architecture is described in Fig. 11, is the use of end-use to power 
demand conversion schemes. This leads to a simple modelling 
way rather than a very precise one complying with the authors’ 
wish that is a satisfactory load curve representation at a 1h 
resolution. 

Concretely, for each domestic activity that implies an 
electricity and/or a water consumption the authors assign a 
standard end-use pattern; each of them is described by a 
series of parameters. These were set so that a realistic average 
consumption is got when they are applied to the simulated 
households. Nominal power and use duration were estimated, 
thanks to a campaign of tests. Activities whose electricity and/ 
or water needs were insignificant have been neglected. This 
modelling approach constitutes an essential feature of Widén 
et al.’s model. 

Thus, five possible schemes have been built up: 


1. either power demand is disconnected from the activity. That is 
the case of cold appliances whose modelling is a constant 
power demand in functioning mode; 

2. or power demand is constant during the activity. This scheme is 
applied to a lot of domestic activities from cooking to ironing 
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and to TV watching. Thus it constitutes a strong assumption in 
the sense that it is applied to very different appliances; 

3. or power demand is constant after the activity. Dishwasher, 
washing-machine and tumble-dryer are concerned by this 
scheme. In fact, the corresponding domestic activity is the 
filling up of the washing appliance. Although power demand 
effectively is required immediately after the end of the 
domestic activity, the authors admit that they widely simplify 
the reality; 

4. or power demand is constant during the activity with addition 
of a temporal constraint. This scheme is applied to the 
majority of activities that imply domestic hot water especially 
when drawings-off occur in the first part of the task. The 
constraint restrains the reuse of already drawn hot water for 
the following activities; 

5. or finally the domestic activity requires a fluctuating power 
demand. The modelling of lighting follows this assumption. 
More precisely, the power demand fluctuates according to the 
available daylight: if somebody is at home in an aware state 
(i. e. not sleeping) the maximum power demand level Pmax is 
applied when the daylight availability L is low (under a Lmin 
level) and the other way around (Pmin When L > Lmax). 


Because of the lack of information concerning the use of DHW, 
the authors had to make use of empirical rules. 

Widén et al.’s model was developed in MATLAB. Load curve 
time resolution can be selected by the program user (from 5 min 
to 1h). Distinction between detached house and flat on the one 
hand and between weekday and weekend day on the other is 
made by the simulation tool. 

As validation, the authors simulate the power demand require- 
ments of the households described in the 1996’s statistical survey 
(that contains among other things time of use data and equip- 
ments characteristics) and compare them with the results coming 
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Fig. 11. Bloc diagram of Widén et al.'s model [26]. 
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from the latest (2006) available load curve measuring campaign. 
Thus, two assumptions are conceded: domestic appliances and 
TOU data have not changed significantly between 1996 and 2007 
enabling the possible use of the input dataset. 

To some extent, the results (individual households’ total load 
curves, i. e. all simulated end-uses aggregated) are in accordance 
with the comparison references notably in terms of load curve 
shape and load peak prediction (high consumption days and 
households). 

Yet, Widén et al. notice that the relevance of the results 
becomes weaker when considering the end-uses separately. For 
instance, the power demand modelling for end-uses that are 
simulated according to the second conversion scheme is more 
satisfactory than these when power requirement and activities 
are disconnected. 

Moreover, the authors carried out other simulations in order to 
make comparisons at an aggregated scale especially on a house- 
holds sample that represents the Swedish population. The model 
shows a good ability to differentiate the average power demand 
shape that characterizes a detached house from an apartment on 
the one hand, a weekday from a weekend day on the other. To go 
further in the comparison, the authors calculated the NVF for 
some end-uses: this shows similarities between measures and 
simulations. Yet when observing more precisely the outputs of 
the simulation tool, some discrepancies can be shown as for 
example an undervaluation of the power demand for a detached 
house during the night. Widén et al. suggest that this is due to the 
miscellaneous end-uses that are not considered by the model. 
Other divergences can be seen as for instance temporal move of 
the modelled peak load for some end-uses compared with 
measurements. 

From TOU data and information about the stock of domestic 
appliances, the authors built up a simulation tool that recon- 
structs load curves for different households’ and days’ typologies. 
With the help of a restricted number of domestic activity to 
power demand conversion schemes, Widén et al. get satisfactory 
results when considering the load curves at aggregated scales 
(wide simulated sample, 1 h time resolution). 

However, the tool is directly modelled on the available time 
of use data, that is to say that it represents the 1996’s situation 
(non scalable modelling). On top of that, these data have been 
collected in diaries that were filled in by the selected house- 
holds’ occupants. The use of the collected data needs a good 
quality of them what is not ensured as underlined by the 
authors. Moreover, there is no modelling of electricity con- 
sumption diversity because it is embodied in the input data. 
As a consequence, the integration of a new end-use seems to be 
difficult because it implies the implementation of the corre- 
sponding time of use scenario in the available diaries. The 
model corresponds to the Swedish housing stock where electric 
heating appears as an exception. That is the reason why this 
end-use (and the thermal issue in general) is not considered in 
the model. 

At the end of [26], Widén et al. list the possible improvements 
of their model: they notably announce the use of a Markov-chain 
and of density functions to generate diversity in terms of human 
behaviour on the one hand, its changes in time on the other. 
These evolutions have been implemented later in Widén and 
Wackelgard’s model. 


3.4.5. Widén and Wéckelgard’s model 

In [30], Widén and Wdackelgard expose the model they devel- 
oped: this is an high time and spatial resolution tool that 
generates domestic activities scenarii. These are then converted 
into electric load curve. 


The main improvement of this model compared with Widén 
et al.’s is the use of a non homogeneous Markov-chain that 
generates diverse domestic activities sequences with high resolu- 
tion levels (1 min time steps, modelling of each household’s 
member). Only the activities that imply electric power demand 
are taken into account by the model: heating and domestic hot 
water supply’? are notably excluded. 

Load curve synthesis follows the same two-step process as in 
Widén et al.’s model: 


1. domestic activity patterns are allocated to each dwelling’s 
occupant; 

2. activities are converted into power demand thanks to conver- 
sion schemes. 


The functional architecture of the model is schematically 
represented in Fig. 12. 

As expected, input data of this model is the same as the 
previous model's: time of use for domestic activities coming from 
statistical surveys and weather data. From these sources of 
information, a Markov-chain (firstly applied to the lighting and 
fully detailed in [31]) has been built up: this is the main feature of 
Widén and Wackelgard’s modelling approach. The model assumes 
that at any simulated time t, each occupant of a household must 
be in one of the three following availability states: 


e absent: the person is not at home; 

e present and inactive: the occupant is at home and sleeps; 

e present and active: the dwelling’s member is at home and can 
carry out each kind of domestic activity. 


Nine domestic activities/end-uses'? are dependent on 
these availability states especially the “present and active” 
that makes them possible: in the Markov-chain, they are repre- 
sented by a code (number). Domestic cold and lighting are 
independent of last availability state and are not driven by the 
Markov-chain. 

Starting with the carrying out of the domestic activity a, at 
time t, the stochastic process calculates the transition probability 
to change to activity daz (a, = dz or a; # az) at t+ 1 with help of the 
input data. Transition probabilities have been calculated at a 1h 
time resolution. Even if load curves can be computed at a better 
resolution, activities’ transition probabilities stay constant over 
one simulated hour. However the probabilities are not fixed along 
the day (non homogeneous Markov-chain), so the daily fluctua- 
tions of the activities are taken into account by the process. 
At each time step, these calculated figures are compared with a 
random number generated by a [0;1] uniform density function. 
According to the sign of the difference, activities are launched 
or not. 

The detailed calculation process is described in [30]. For few 
cases, the authors adapted it to ensure its coherence (especially 
for night periods). 

As in Widén et al.’s model, domestic activity to load curve 
conversion schemes have been built up with help of measured 
load curves. Moreover these have been refined compared with 
those in [26] especially in three points: 


1. appliances’ power demand can occur before, during, just after 
the corresponding domestic activity or even totally discon- 
nected from it; 


12 Very little electricity demand is dedicated to these end-uses in the Swedish 
housing stock. 

13 Away, sleeping, cooking, dish-washing, washing, TV, computer, audio 
appliances and other devices. 
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Fig. 12. Bloc diagram of Widén and Wackelgard’s model [30]. 


2. unitary power demand pattern and stand-by power level are 
characteristic of each appliance; 

3. appliance sharing between household’s occupants is taken into 
consideration. 


Concerning this last point, the authors underline that this 
kind of data are not recorded in the time of use diaries. 
Thus a supplementary assumption is made by Widén and 
Wackelgard: a few appliances are totally shared by all dwelling’s 
members (e.g. the fridge) whereas some other are systematically 
used by a single occupant. According to the device status, 
the related domestic activity can be carried out by several 
people at the same time (leading to a single power demand 
cycle) or not (the activity can only be conducted by a single 
person and it systematically leads to a power demand 
pattern). 

In [30], more details are given about the modelling of each 
activity related appliances. In fact, assumptions should have be 
done in order to adapt the modelling method to each domestic 
device (respectively each appliance related to the same domestic 
activity). 

So as to validate their model, Widén and Wackelgard pro- 
ceeded with two kinds of modelling tasks and comparisons taken 
the Sweden Energy Agency’s data as references: first they tested 
the model on a small scale housing stock (14 dwellings) in order 
to evaluate the ability of the model to generate diversity. Then 
they compared the generated load curves at a more aggregated 
scale (two hundreds simulations for a similar dwelling type 
during a week - five weekdays and two weekend days, two 
dwelling types considered detached house/flat) so as to know 
if the model provides a good approximation of the energy 
consumption with regard to the statistical national energy 
consumption data. 


To compare simulated and measured load curves, the authors 
chose the two following criteria: 


e the load factor Fn, Pmax related to the maximal power demand 
for a set of n clients during an interval At: 


Pmes,n(At) 


Fn,Pmax( At) = Prmaxn(Ab) 
max,n 


(10) 
with Pmesn and Pmaxn respectively the mean and maximal 
power demand required by the n consumers over the time 
period At 

e the diversity factor Ky, related to n electricity users and the 
time slot At: 


Xj- 1 Pmaxj(At) 


Krna = Prax OO 


(11) 
where Pmaxj is the maximal individual power demand value 
for the client j and Pmax,n is the maximum demand of the n 
consumers considered in the whole. 


Both equations give a summarized analysis when considering 
different power demand profiles and thus explain their frequent 
use in such a study. 

In the whole, comparisons provide the same kind of conclu- 
sions that have been made for Widén et al.’s model: the simula- 
tion tool gives satisfactory estimations of domestic load curves, 
generates sufficient differences in terms of timing and magnitude 
of the power demand between weekdays and weekend days on 
the one hand, between house and flat on the other hand. 

However some improvements can be carried out such as a 
better modelling of some end-uses, a refinement in terms of 
domestic activities and a more suitable distribution of the 
domestic equipments according to the dwellings’ properties. 
Carrying out these evolutions necessarily requires more detailed 
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input data. These come from surveys and their precision directly 
depends on the attention of the selected people sample. More- 
over, the model is only relevant for short term forecasting since 
time of use data always represents the past situation. This is 
notably what causes the little ability of the model to simulate the 
Information and Communication Technologies (ICT) load curve 
(the computer is the best example). Thus to maintain a satisfac- 
tory precision of the model, updated time of use data have to feed 
the simulation tool imposing regular to continuous surveys. 
Although heating is not considered in this model, the authors 
underline that its integration is feasible. Yet, without apparent 
help of a building simulation tool, we can wonder about the 
precision of this end-use modelling. 


3.4.6. Richardson et al.’s model 

In order to quantify the impact of “low carbon measures”, such 
as demand response and CHP, on a local electricity distribution 
network, Richardson et al. developed a model presented in [32]. 

The authors adopted a “deep” bottom-up approach since the 
simulated elements are the specified household’s domestic equip- 
ments and the households’ members whose individual behaviour 
modelling is based on activity diaries built up from surveys. Thus 
Richardson et al. are able to reconstruct the 1 min electric load 
curve corresponding to a few number of dwellings located in the 
same area. According to the authors, this simulation time resolu- 
tion has been considered as the best compromise between the 
amount of produced results and the required precision level at a 
small spatial scale. 

The structure of the simulation tool is schematically drawn in 
Fig. 13. 

As Fig. 13 underlines, Richardson et al. made use of the results 
of a time of use survey in order to get availability scenarii 
(presence at home and awake state). According to the authors, 
the occupant’s state “at home and active” is the most interesting 
because in this case electric devices can be run. In the model, this 


kind of information is coded by a binary variable which equals 
1 when an occupant is at home and active, 0 otherwise. An 
aggregated availability profile is obtained when adding each 
occupant’s availability variable. 

Yet the active occupation only corresponds to the necessary 
condition to use an equipment, the next step to construct a load 
curve is to know the kind of device that is run, the use duration 
and the start times. That is the reason why Richardson et al. have 
made use of daily activity profiles (probability density over the 
day at the diaries’ time scale) elaborated from the TOU survey for 
seven different domestic activities. Here we recognize the same 
input material as in Widén and Wéackelgard’s method [30]. 
However these elements differ from the Swedish model because 
they are “static”: in Richardson et al.’s model, daily activity 
profiles are only functions of the day-type (weekday/weekend 
day) and the number of occupants in the simulated household 
(from 1 to 5). Using such kind of data enables a realistic human 
behaviour modelling on condition that first diaries have been 
written down carefully and secondly the people sample is large 
enough to “smooth” domestic habits. 

So as to provide a satisfactory device modelling, Richardson 
et al. used statistical information (penetration rates, devices’ and 
households’ annual consumption) and some measured power 
demand patterns - all input data are available in [33]. When this 
last type of information was unavailable for the appliances, 
constant power demand has been assumed. They sometimes 
had to adjust them in order to reach the average East Midlands’s 
household that constitutes the target annual domestic 
consumption. 

Domestic activities can imply the need of electricity without 
forcing it. That is the reason why Richardson et al. had to make 
connections between activities and appliances: according to the 
authors’ classification of the considered domestic activities 
(watching TV, cooking, laundry, washing/clothing, ironing, clean- 
ing, other activities), some links are obvious because one activity 
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Fig. 13. Bloc diagram of Richardson et al.’s model [32]. 
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corresponds to a single device (watching TV «television, 1-1 
relation), other are more “ambiguous” (cleaning < ?, 1-n link). 
However, the last case does not mean that when carrying out an 
activity linked with n electric devices, all n appliances are turned 
on as long as the activity is conducted. It only allows a possible 
simultaneous use of the same activity related appliances. 

“Other activities” related appliances have been put together by 
the authors because: 


e either any of the other activity profiles did not match with 
them (e.g. telephone). Thus they are only dependent on the 
human presence at home, 

e or their functioning is continuous (e.g. cold appliances) and 
their modelling is not function of the availability at home, 

e or finally their use is much more linked with an exogenous 
parameter such as the outer temperature what is the case of 
heating. 


When simulating a sample of households, their set of domestic 
appliances has been randomly assigned with respect of the 
national penetration rates (multi-equipment is taken into 
account). During a simulation, an electric device can be in two 
different states: turned on and turned off. For some appliances, 
last state imply a standby consumption which is taken into 
consideration in the load curve modelling. Moreover, in Richard- 
son et al.'s opinion, sharing of equipment and use of linked 
devices (e.g. washing-machine/tumble-dryer) can be “read” directly 
from the time of use data that have been arranged in the adequate 
form to serve as input of the model. 

To calculate a load curve at each time point, four steps are 
required: 


1. according to the selected appliance, the number of active 
people at the simulated time step and the type of the day, 
the corresponding daily activity profile is chosen; 

2. the calculated time point value is read on the activity profile; 

3. then the device starting probability has to be calculated. It is 
done by multiplying the value obtained in step 2 with a scalar 
factor. In fact, this last figure ensures the energy accuracy of 
the model because it relies on the mean annual number of uses 
for each equipment (device energy consumption/mean energy 
consumption per cycle); 

4. finally, the device launching probability is compared to a 
random number generated thanks to a [0; 1] uniform density 
function. If the probably is greater than the random number, 
the appliance is launched. 


In order to validate the results of their model, Richardson et al. 
conducted an in-depth multi-criteria study. Comparisons have 
been carried out considering following aspects: 


è annual and mean daily energy consumption; 
e daily load curves; 

e 1 min power demand volatility; 

e diversity of the power demand; 

e load duration curve; 

e load factor. 


As reference data, Richardson et al. make use of 1 min load 
curves measured in 22 dwellings. These data are totally indepen- 
dent of the values that serve as model inputs. 

As reference data, Richardson et al. make use of 1 min load 
curves measured in a sample of 22 dwellings located in the East 
Midlands area. These validation data are totally independent of 
the values that serve as model inputs. 


The followings points sum up the main conclusions made by 
the authors: 


e electric heating has been excluded from the validation study 
because of a too low penetration rate; 

è one year simulation for 22 households shows that calculated 
energy consumption matched well with the validation data in 
average but its standard deviation was less than the real one. 
Assumptions of the model and input data biases can be the 
causes; 

e in terms of average daily consumption, inter-month variations 
are correctly modelled but seasonality influence is under- 
evaluated. Lack of precision in the diaries and non modelling 
of the heating auxiliaries can explain this observation; 

e when considering the mean daily calculated load curve, the 
model under-estimates power demand during the night and 
the morning peak occurs later than in the reality. The fact that 
the measured households sample does not match the national 
case can be responsible for this difference; 

e the model shows a little ability to represent little and 
large power demand transitions.'4 To the authors’ opinion, 
modelling simplifications are largely responsible for these 
observations; 

e diversity factor and After Diversity Maximum Demand (ADMD) 
values from the model are close to the real indicators; 

e load duration curves from reality and model are quite similar. 


As a conclusion, Richardson et al.’s model provides high spatial 
and time resolution load curves thanks to detailed time of use 
data. Its strengths and weaknesses have been identified by the 
authors and mostly listed in the previous paragraphs. We would 
underline the unknown concerning the modelling of heating. In 
[32], Richardson et al. do not mention the use of a building 
simulation tool which seems to us necessary to build up a 
relevant load curve at a fine time step. Concerning the appliances 
modelling and more precisely the nominal power values of the 
equipments, only mean values have been used in the model. A 
statistical distribution can generate more variations and diversity 
between the simulated households. Moreover, the model is 
restricted to short-term forecasting because new end-uses cannot 
be integrated. Finally, households’ socio-economic characteristics, 
price of energy and DSM measures are load curves influences that 
are not included in Richardson et al.’s model. 


4. Statistical-engineering model: Train et al.’s model 


Previous models either follow an upward (bottom-up) or 
downward (top-down) process that both present advantages 
and limits. From this assessment, the main motivation to develop 
an other kind of model is to put to use all benefits of each 
approach. That is what Train et al. [34] tried to do in building up a 
Statistical Adjusted Engineering (SAE) model. 

The methodology corrects the end-use load curves simulated 
with an engineering approach (so outside the model) thanks to 
statistical coefficients obtained through measured load profiles. In 
terms of input data, Train et al. made notably use of a PG&E 
metering campaign [35] that provided them total and by end-use 
load curves for about 800 domestic clients. During this survey 
information concerning the households (income, dwelling size 
area, number of occupants, appliances penetration rates, etc.) has 


14 Transition = |P(n+1)—P(n)| where P is the power demand value at the time 
step n. 


A. Grandjean et al. / Renewable and Sustainable Energy Reviews 16 (2012) 6539-6565 6559 


Operational elements 
of the model 


Output data 


Geometric and thermal 


dwelling characteristics 


Weather data 
Appliances penetration rates 


Presence at home scenarii and 
Time Of Use data 


Socioeconomic data: 
household income and number 
of occupants 


Total load curves — 
800 clients 


End-use load curves — 


2 monitored appliances 
in 800 clients 


Heating load curve for 
each dwelling (60 min) 


Engineering load curve for 
each appliance (60 min) 


SAE end-use load curves 
(60 min) 


End-uses aggregation 


SAE total load curve 
(60 min) 


Fig. 14. Bloc diagram of Train et al.’s model [34]. 


been collected. Moreover engineering loads have been calculated 
by ADM Associates using the method described in [36]. 

Thus engineering power demand values correspond to the 
explanatory variables of a conditional demand model into which 
they are put. The general functioning of Train et al.’s model is 
shown under a bloc diagram in Fig. 14. Results are given at the 1h 
time resolution. 

As it is said in [35], statistical adjustments on load curves can 
occur at each time step (1 h), for each end-use i and/or for each 
customer j. However it does not seem to be relevant to adjust 
engineering load curves on these three sources of influence 
because it would totally “hide” the simulated part of the load 
curve results. 

That is the reason why the authors constructed four SAE 
models with progressive adjustments. These are applied to one 
or several time periods | (I¢[1;1]) defined within the day on the 
one hand and to take into account the customers’ characteristics 
on the other hand. 

In the first defined SAE method, engineering load curves are 
adjusted with a correction factor zeta; that is specific for each 
end-use i (zeta; = zeta;(i)). Yet this coefficient is constant over the 
day (only one period is defined) and it does not take into account 
the customers’ characteristics. In the second SAE model, these 
characteristics are taken into consideration for the adjustment of 
the load curves (zeta; = zeta;(i,j)). However, once again only one 
period is defined for the adjustment that is to say that zeta; is not 
time-dependent. In the third SAE approach, four periods | are 
defined and the adjustments are different according to the end- 
use and the period but without considering the customers’ 
specificities (zeta; = zeta;(i,l)). Finally the last SAE model variant 
is the combination of model 2 and 3: engineering load curves are 
corrected according to the end-use i, the customer j and the 
period l (zeta; = zeta;(i,j,l)). 


Whatever the adjustment level, Train et al.’s model is looking 
for estimating the total power demand of customer j at time t, 
P(t), in rectifying the simulated power demand value at t for the 
end-use i of the customer j, R; j6,” thanks to the correction factor 
zeta;;jı and an error u,(t) according to Eq. (12). 


n L 
P= X $ zeta;iji - Rig(O) - Dij+ u0) (12) 


i=1l=1 


In few words, the authors observe that different adjustment 
periods provide a better correction of the load curves calculated 
with engineering models. However, introducing demographic 
variables (customers’ characteristics) in the model do not seem 
to be of great interest when considering the gain in terms of 
results accuracy. 

A wider evaluation and validation work has been conducted by 
Train in [35]. This publication mentions an extended version of 
the fourth model (adjustment on nine periods, differentiation in 
terms of type of day, adaptation of the model according to the 
season summer/winter, etc.) but conceptually it is the same 
method than this explained above. 

In the whole, the assessment of the model has lead to general 
conclusions that we previously mentioned as for instance the 
limits of the engineering models to take into account the time of 
use and the regulation of thermal appliances. On top of that, Train 
compared the SAE, engineering and measured daily load curves 


13 n end-uses are considered - a dummy variable Dj; indicates the presence or 


the absence of each of them in the customers’ homes. 
16 That depends on the end-use i and/or the consumer j and/or the adjustment 
period | according to the SAE model variant. 


Table 2 


Reading grid (first part) of the studied load curve reconstitution models. 


Type of model 


Authors [Publication] 
Purpose of the model 


Inputs of the model 


Outputs of the model 


Modelled 
end-uses 


Nature of the generated 
behavioural diversity 


Model validation 


Deterministic statistical 


disaggregation 
(top-down) models 


Statistical random 
bottom-up model 


Time of use based 
bottom-up models 


Aigner et al. [10] 


To get end-uses load curves 
without monitoring 
Bartels et al. [12] 


Simulation of the impact of various 
scenarii on the load curve at the 
regional level 

Decision support for the planning 
of the capacities of the utilities 


Yao and Steemers [14] 


Support for the design of energy 
systems including renewables 
Support for the prediction of the 
load curve for a selected community 


Walker and Pokoski [23] 


Prediction of the load curve as 
support for the planification of new 
power generation capacities 
Capasso et al. [25] 


Support for the utilities for the 
evaluation of power management 
measures and DSM policies 


Armstrong et al. [29] 


To get electricity consumption 
profiles for the modelling of 
micro-cogeneration devices 


Assumptions concerning the non 
functioning of a series of equipments 
at certain hours within the day 
Measured daily load curves at the 
household level 

Evolution scenarii concerning 
technological aspect, socio- 
economical aspect, appliances 
penetration rates, number of 
customers 

Measured daily load curves at the 
household level 


Occupation scenarii 


Statistical daily consumptions by 
end-use and appliance 


Occupation, activities and 
end-use scenarii 


Unitary load cycles (from modelling 
and apparently from measures) 


Occupation and activities scenarii 


Constructed unitary load cycles 


Consumption available at the 
household level (with an ambiguous 
temporal restitution of these data) 
Statistical mean annual consumption 
per appliance (from measures but 
with an ambiguous temporal 
restitution of these data) 

Occupation and activities scenarii 


Measured and constructed unitary 
load cycles 


Mean annual consumption available 
at the household level 

Mean annual consumption of the 
appliances 


Daily load curves returned at a 60 min time step for nine 
domestic appliances at the household level 


Daily load curves per appliance returned at a 60 min time 
step for the regional residential class considered. The load 
curves are calculated according to a scenario (worked day 
or not worked day for a selected month) 


Daily multi end-uses load curves returned at a 1, 5, 15 or 
30 min time step corresponding to, from a single 
household, to an entire community 


Daily load curves returned at a 15 min time step for a 
household (or a sample of households) that is specified 
with taking into account of the psychological and 
behavioural influences 


Residential daily load curve returned at a 15 min time step 
for a geographical area thanks to the synthesis of multi 
end-uses load curves at the household level 


Daily load curve for electricity specific appliances returned 
at a 5 min step for typical households 


Specific electricity 
appliances, heating, 
cooling, DHW 


Specific electricity 
appliances, heating, 
cooling, DHW 


Specific electricity 
appliances, heating, 
DHW 


Specific electricity 
appliances, DHW 


Specific electricity 
appliances 


Specific electricity 
appliances 


Deterministic 


Deterministic 


Mathematical random 


Scenario based 
probability 


Scenario based 
probability 


Scenario based 
probability 


No validation 


On total load 
curves 


On total load 
curves 


On total load 
curves 


On total load 
curves 


On reconstructed 
load curves for 
specific electricity 
appliances 
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Table 3 


Reading grid (second part) of the studied load curve reconstitution models. 


Type of model 


Authors [Publication] 
Purpose of the model 


Inputs of the model 


Outputs of the model 


Modelled 
end-uses 


Nature of the generated 
behavioural diversity 


Model validation 


Time of use based 
bottom-up models 


Probabilistic 
empirical 
bottom-up models 


Widen et al. [26] 


Determination of the households’ 


energy consumption 
Evaluation of future changes 


(modifications of the behaviour, 
energy efficiency improvement) 
that impact on the electricity use 


Widén and Wackelgard [30] 


Support for studies on electricity 
production at the household level 


Richardson et al. [32] 


To anticipate the impact of “low 
carbon” measures on the load 


curve of a local distribution 
network 


Stokes [17] 


Evaluation of the impact of the 


decentralized electricity 


production in the case of a low 


voltage network 


Paatero and Lund [21] 


To get accurate domestic 


electricity consumption data to 


identify the impact of DSM 
measures 


Occupation and activities scenarii 


Constructed unitary load cycles 


Occupation and activities scenarii 


Measured daily load curves per 
appliance 
Constructed unitary load cycles 


Occupation and activities scenarii 


Measured unitary load cycles 


Statistical yearly consumptions 
per appliance based on studies 


Daily load curves measured at the 
household level 


Daily load curves measured at the 
appliance level 


Unitary load cycles constructed 
from measures 

Yearly consumptions per end-use 
coming from surveys 

Daily load curves measured at the 
household level 


Unitary load cycles constructed 


from measures 


Daily consumptions measured at 
the household level 


Daily load curves returned at a 5-60 min time step 
and DHW draw profiles (ambiguous temporality) at 
the occupant level. Aggregation of previous results to 
get those corresponding to a household or a little 
community 


Daily load curves returned at a 1-60 min time step 
for a household and possible aggregation to get 
results at a larger scale 


Daily load curves returned at a 1 min time step for a 
household and possible aggregation for a set of 
households 


Daily load curves returned at a 1 min and 30 min 
time step for respectively a specified and “mean” 
household and possible aggregation at a community 
scale 


Daily load curve at a 60 min time step for a household 
and possible aggregation at a larger scale (thousands 
of households) 


Specific electricity 
appliances, DHW 


Specific electricity 
appliances 


Specific electricity 
appliances, 
heating, DHW 


Specific electricity 
appliances, 
heating, DHW 


Specific electricity 
appliances 


Scenario based 
determinism 


Scenario based 
probability 


Scenario based 
probability 


Empirical probability 


Empirical probability 


On load curves per 
end-use 


On DHW draw profiles 


On total load curves 


On daily electricity 
consumption per end-use 
On yearly consumption for 
DHW 

On total load curves 


On the household’s yearly 
consumption 

On yearly consumption per 
end-use 

On total load curves 


On total load curves 


On total load curves 


On household's yearly 
consumption 
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Table 3 (continued ) 


Nature of the generated Model validation 


behavioural diversity 


Modelled 


Inputs of the model Outputs of the model 


Authors [Publication] 


Type of model 


end-uses 


Purpose of the model 


Consumptions per appliance in 
proportion with the dwelling 


consumption (data from studies 


without a precise temporality) 


On total load curves 


Daily load curves returned at a 60 min time step per Specific electricity Statistical 


appliance and for a specified household (or a group of appliances, 


households) 


Measured daily load curves 


Train et al. [34] 


Hybrid model 
(statistical- 


available at the household level 


heating, cooling, 


DHW 


engineering) 
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Appliances daily load curves from 


modelling 


To take advantage of the 


statistical and the engineering 
models to obtain end-uses load 


curves 


with help of the statistic S defined in Eq. (13). 


VEE (LO-MO? 
S= a MD -24 (13) 
where L(t) and M(t) are the 1 h power demand value averaged on 
the sample of simulated customers that is obtained by the SAE (or 
the engineering) method and measured respectively. In fact, this 
is the square root of the NVF that we previously have met. 

More precisely, Train focused his assessment on the compar- 
ison between SAE and engineering loads in order to quantify the 
impact (improvement and/or deterioration according to the 
considered end-use) of the statistical adjustment. The conclusion 
of the author is that the SAE model improves the engineering 
results if and only if physical models do not integrate empirical 
knowledge (so from measures) or the structure and/or the size of 
the sample used for the adjustment is (are) not optimal in terms 
of representativeness. 

To conclude, Train et al. developed an hybrid model that takes 
advantage of the traditional load curve modelling techniques. 
Measured load curves calibrate power demand profiles obtained 
by engineering models. On top of that, this method allows to get 
round the diversity modelling because it is contained in the 
adjustment coefficients. 

However this approach is only relevant for short term fore- 
casting because it reproduces past power demand profiles what is 
the case of some other models. Changes and evolutions in the 
housing stock cannot be implemented and simulated in this 
model. On top of that, the miscellaneous category encompasses 
a lot of equipments. The corresponding load curve can be well 
estimated in the whole but what about each of the appliance in 
this field? 


5. Cross analysis 


In this section, the authors sum up and compare the models of 
the literature that were described in the previous paragraphs. 
First of all, a general overview of them is given under a table 
form. This contains the in-our-opinion most important criteria 
that characterize a load curve reconstitution methodology. 
Then certain comparisons are proceeded with focusing on 
the advantages and the limits of each methodology. Finally, we 
plot the results of this cross analysis through a “pseudo 3D” 
representation. 


5.1. Reading grid 


To construct this reading grid, we began to list the to-our- 
mind most relevant characteristics for defining a load curve 
reconstitution model: these criteria correspond to the columns 
of the table and the models are arranged in rows. 

Then we completed each cell of the table thanks to a restricted 
number of modes so as to facilitate the reading. This way we 
obtained Tables 2 and 3. 


5.2. Focused comparisons 


5.2.1. Modelling the diversity 

As we previously said, estimating the electric power demand 
for the domestic sector implies the modelling of the diversity. The 
above descriptions of the methodologies show that various 
strategies have been developed by the authors to represent it. 
Some of the models only “catch” the diversity directly from the 
input data without clarifying it. Other models generate it totally 
or partially thanks to more or less elaborated random processes 
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(probability distribution functions, Monte Carlo extractions, Mar- 
kov chains, etc.) or statistical approaches (correction factors). 

When only considering the models that represent the diver- 
sity, Walker and Pokoski’s work is the to-our-knowledge first one 
which separately takes into account the presence at home of 
the households’ members on the one hand (necessary but not 
sufficient condition to reconstruct the domestic load curve) and 
the proclivity for carrying out certain domestic activities (or using 
appliance(s)) on the other. On top of that and like the models of 
Capasso et al. [25], Widén et al. [26], Widén & Wackelgard [30] 
and Richardson et al. [32], each household’s member is characterized 
such a way. This constitutes a real modelling refinement. 

However, in order to feed models that follow such an 
approach, precise information have to be available. Among other 
data, the authors made use of time of use surveys which are only 
representative of a selected households sample at a certain time. 
In order to ensure the representativeness of the models using TOU 
data, these must be periodically updated which supposes perma- 
nent and so expensive studies. 

In some models, a consumer mean behaviour has been 
defined and applied to the simulated number of dwellings. Here 
is open the question of the likelihood of this construction 
because a mean consumer only has sense with regard to one 
of its characteristics or some of them. As the behaviour is 
influenced by too many parameters, defining a realistic mean 
one does not seem relevant. 

For energy purposes, building up a mean consumer is viable 
because it is possible to calculate the mean consumption of a 
consumer sample over a certain duration: the mean consumer is 
the one who uses the calculated average energy consumption. 
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In terms of power demand, a mean consumer has no sense; 
because of the diversity in a customers sample, that among 
other things comes from behavioural and technical influences, 
there is every chance that nobody has the same load curve as 
the After Diversity Mean Load Curve (ADMLC). It is calculated 
by addition of each individual load curve and then division by 
the number of aggregated consumers. However what can be 
defined with regard to the considered sample is the marginal 
consumer that is to say the consumer that will be added to the 
original sample. Its most likely representative load curve would 
be the ADMLC. 

Another possibility to model the diversity through the human 
behaviour is to calculate from input data transition probabilities 
from a domestic activity to another one, and this, for all activities. 
When integrating this within a Markov Chain, a random influence 
completes the variety that is already contained in the input data. 
However the link between the activities and the use of appliances 
remains to be set up what is not trivial. 

To our opinion, it is possible to set a framework for the 
modelling of the diversity that can be less dependent on the input 
data. With help of a sufficient knowledge of the simulated housing 
stock, logic rules and constraints can be defined and applied in 
order to generate diversity. 


5.2.2. Allocation of appliances functioning cycles within a simulated 
day 

Different approaches exist in order to construct scenarii that 
define the use of appliances at the daily level. 
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Fig. 15. Cross analysis diagram of the load curve models. 


6564 A. Grandjean et al. / Renewable and Sustainable Energy Reviews 16 (2012) 6539-6565 


The first one is to calculate an end-use average daily or weekly 
consumption from annual consumption figures and to launch a 
corresponding number of appliances cycles according to the 
unitary consumption of each end-use. However, this method 
can lead to an unrealistic daily usage of each appliance as for 
instance a ¢ washing-machine cycle to launch each day. 

The second one is to make use of the daily end-uses (or domestic 
activities) probability distribution functions constructed thanks to 
time of use data. In order to reach the energy consumption target, 
the number of end-uses starts must be controlled. This way is 
limited in the sense that in the better case it reproduces one 
measured reality and is representative of it. Changes in the 
end-use behaviour can not be simulated. 

To our mind, we think that a method based on the knowledge 
of appliances use frequencies at different time scales can be 
more realistic than the previous ones. It will rely on the 
definition of types of day and households. When considering 
the week as the modelling framework, the sequential distribu- 
tion of the launching of end-uses cycles for each day will 
take into account these influences but also the household’s 
availability at home duration and the number of end-uses cycles 
already started in the considered week. On top of that, it would 
be necessary to constitute some links between appliances in 
order to improve the realism of simulations. For instance, the 
use of tumble-dryer is connected to the use of washing-machine. 
Moreover, it will be required to construct an end-uses ranking 
for each simulated day so as to launch them according to a 
priority order. 


5.3. Positioning the models on a pseudo 3D plot 


As we said in the introduction, we decided to represent on 
a chart the positioning of the models. In order to facilitate 
the reading of the representation, we restricted the number of 
regarded criteria: thus we only considered the range of appli- 
cation of the model (so the simulated appliances), the way 
the model (re)constructs the diversity and finally its time 
resolution. 

In Fig. 15, the models are represented by bubbles whose 
diameters are inversely proportional to their time resolution. 
That is to say that a big bubble corresponds to a model that 
provides load curves with a fine time accuracy. 

In the north-east corner of Fig. 15, a lack of models is notice- 
able: to our knowledge, there is no model that simultaneously 
provides load curves for all domestic appliances (inclusive ther- 
mal ones), at a fine time scale and with an advanced methodology 
to represent the diversity. 


6. Conclusion 


In this paper, we review 12 domestic load curve models. 
In the first part, we adopted a systematic analysis of each of 
them with focus on the finality of the model, its inputs/ 
outputs, the modelled end-uses, the way diversity is taken 
into account and finally what kind of validation has been 
carried out. On top of that we built-up a synthetic block- 
diagram of each approach. In the second part of the paper, we 
considered them together carrying out a transverse analysis. 
A reading grid and a pseudo 3D plot was established at the end 
of the review. 

This work has been voluntarily limited to a little number of 
models because of the willingness to analyse each of them in 
details. Other existing techniques for the load curve modelling, 
such as neural networks or agent-based method, have not been 
dealt with in this review. Further research has to be carried out to 


list, evaluate and analyse all possible methodologies for the 
reconstitution of the residential load curve. This purpose seems 
to us of great importance because of the undergoing radical 
transformation of the dwellings, domestic end-uses and human 
habits and behaviour. 
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